PSR J0045-7319: A DUAL-LINE BINARY RADIO PULSAR 
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ABSTRACT 

Binary radio pulsars are superb tools for mapping binary orbits, because of the 
precision of the pulse timing method (e.g Taylor & Weisberg 1989). To date, all orbital 
parameters for binary pulsars have been derived from observations of the pulsar alone. 
We present the hrst observations of the radial velocity variations due to the binary 
motion of a companion to a radio pulsar. Our results demonstrate that the companion 
to the Small Magellanic Cloud pulsar PSR J0045-7319 is the B1V star identified by 
Kaspi et al. (1994). The mass ratio of the system is 6.3 ± 1.2, which, for a neutron star 
mass of 1.4 M©, implies a mass of 8.8 ±1.8 M & for the companion, consistent with the 
mass expected for a B1V star. The inclination angle for the binary system is therefore 
44° ±5°, and the projected rotational velocity of the companion is 113 ± 10 km s -1 . The 
heliocentric radial velocity of the binary system is consistent with that of other stars and 
gas in the same region of the Small Magellanic Cloud. 

Subject headings: Pulsars:individual PSR J0045-7319, Stars: fundamental parameters, 
Stars: binaries: spectroscopic 
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The only known pulsar in the Small Magellanic 
Cloud, PSR J0045— 7319, was discovered in a system- 
atic search for radio pulsars in the Magellanic Clouds 
(McConnell et al. 1991). It was recently shown to 
be in a 51-day binary orbit of eccentricity 0.8 (Kaspi 
et al. 1994). The pulsar's mass function of 2.17 Mq 
implies a minimum mass for the companion of 4.0 Mq 
for a 1.4-Mq pulsar. Since the maximum mass of a 
neutron star is ~ 3.0 Mq (Friedman, Ipser & Parker 
1984, Friedman et al. 1988), the companion must be 
either a black hole or a massive non-degenerate star. 
The B1V star identified at the position of the pulsar 
is a strong candidate for the companion (Kaspi et al. 

1994) . 

For a B1V star companion, several observable ef- 
fects are expected, since the pulsar approaches to 
within approximately 5 stellar radii from the star at 
periastron. At radio wavelengths, effects that would 
vary with orbital phase include dispersion, scattering 
and absorption of the pulsed emission. Not one of 
these effects has been observed, however systematic 
frequency-independent timing residuals with respect 
to a Keplerian orbit have been detected (Kaspi et al. 

1995) . Optically, the star should show radial veloc- 
ity variations with an amplitude of approximately 30 
km s _1 . It has also been suggested that the compan- 
ion may exhibit flux variations due to stellar pulsa- 
tions excited resonantly by the eccentric orbit (Ku- 
mar, Ao & Quataert 1995). While the 0.01 magni- 
tude amplitude of this effect will be very difficult to 
detect, such small variations have been detected in 
other objects (Balona 1995). In contrast, if the true 
companion is a black hole and the B1V star is not as- 
sociated, none of the above effects should be observed 
(Lipunov, Postnov & Prokhorov 1995). 

To determine the nature of the companion to 
PSR J0045-7319 observations of the B1V star were 
conducted to search for radial velocity variations. 
Spectra at a resolution of 0.55A per pixel, cover- 
ing 500A centered on 3900A, were obtained at nine 
epochs using the Australian National University's 2.3- 
m telescope at Siding Spring, NSW, Australia. The 
[On] doublet at 3728A from the surrounding Hn re- 
gion was used as a radial velocity standard. Helium- 
argon arcs were used to determine the dispersion. Ra- 
dial velocities were obtained by Fourier cross correlat- 
ing spectra from each epoch with a spectrum obtained 
by summing the spectra from all epochs. The spec- 
tral region used contained the Balmer lines between 
the Balmer jump and H8. 



With the binary period, longitude and epoch of pe- 
riastron, and eccentricity determined from radio ob- 
servations (Kaspi et al. 1995), the systemic radial 
velocity and amplitude of the companion's radial ve- 
locity variation were determined using a least-squares 
fit to the observed velocities. The fit to the radial 
velocity curve is shown in Figure 1 along with the 
residuals and the pulsar's radial velocity curve. This 
fit gives an apparent semi-major axis for the compan- 
ions orbit of 27.7 ± 5.0 lt-sec and hence, a mass ratio 
of 6.3 ± 1.2. Since all measurements of the mass of 
neutron stars are consistent with 1.4 Mq (Thorsett 
et al. 1993), it assumed that the mass of the neu- 
tron star in this system is also 1.4 Mq. Using the 
pulsar mass function of / = 2.17Mq, implies a com- 
panion mass of 8.8 ± 1.8 Mq and an inclination angle 
for the binary of 44° ± 5° . The derived mass is close 
to the typical mass of a B1V star, ~11 Mq (Ander- 
sen 1991). Using the Stefan-Boltzmann law, with an 
effective temperature of 24000 ± 1000 K and lumi- 
nosity of 1.2 x 10 4 Lq, the radius of the companion is 
6.4 ± 0.7 R . 

To test the stability of the fit, the fitting procedure 
was repeated with the point just after orbital phase 
zero removed, resulting in an increase of the mass 
ratio of only 12%. The probability that the results 
shown here are due to random radial velocity varia- 
tions in an isolated B star is very small, considering 
the close positional coincidence (Kaspi et al. 1994), 
the good agreement of the orbital phases of the radio 
and optical radial velocity curves, and the amplitude 
of the optical radial velocity curve which gives a plau- 
sible estimate for the mass of the B star. If we assume 
the points are randomly distributed about a mean 
value, the reduced chi-squared of the prefit data is 
4.4, while that of the postfit data is 0.9. For 8 degrees 
of freedom, the probability of obtaining a chi-squared 
as high or higher than 4.4 is less than 0.001, while 
the probability of obtaining a chi-squared as high or 
higher than 0.9 is 0.5. The accuracy of the measure- 
ments here could be improved with better resolution 
and higher signal to noise ratio. 

A high rotational velocity for the companion might 
be expected as evolutionary histories for this type of 
system involve mass transfer from the pulsar progeni- 
tor to the companion, which serve to spin it up (Bhat- 
tacharya & van den Heuvel 1991). Using the He I 
lines, the projected rotational velocity for the com- 
panion is 113 ± 10 km s _1 . Such velocities however 
are common amongst main-sequence Bl stars. As- 
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compares well with the previously measured value of 
131 ± 3 km s" 1 (Smith & Weedman 1973). The bi- 
nary system's observed heliocentric velocity of 117 ± 
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Fig. 1. — Top: Radial velocity data and fitted curves 
for the radio observations (Kaspi et al. 1995) (large 
amplitude curve) and for the companion, from op- 
tical observations (small amplitude curve). Bottom: 
Residuals of the two-parameter fit to the optical radial 
velocity variations of the companion. In both panels 
the error bars shown are ± one standard deviation. 



suming that the spin and orbital angular momenta 
are aligned the rotational angular velocity of the com- 
panion at its equator is v/r = 3.1 ± 0.3 radians per 
day. When compared with the orbital angular veloc- 
ity at periastron of 2tt_P j " 1 (1 + e)° 5 (1 - e)" 1 5 = 1.96 
radians per day, this indicates that synchronisation is 
unlikely to have occured, although it cannot be ruled 
out given the uncertainty in the B star's inclination 
relative to the orbit. 

The radial velocity of the companion was found to 
be -10 ± 7 km s _1 relative to the surrounding nebula 
NGC 248 (N13A + N13B) (Henize 1956). The helio- 
centric radial velocity of 127 ± 7 km s" 1 for NGC 248 



7 km 



is consistent within the uncertainties of the 



167 ± 40 km s 1 estimate from low dispersion spectra 
(Kaspi et al. 1994). The observed heliocentric veloci- 
ties of both the PSR J0045-7319 binary and NGC 248 
are somewhat less than the mean heliocentric radial 
velocity for the Small Magellanic Cloud (SMC) of 160 
kms -1 (Martin, M aurice & Lequeux 1989). However, 
when one considers the radial velocities of the four 
components of the SMC, -45, -28, 9 and 30 km s _1 
(relative to the mean) (Martin, Maurice & Lequeux 
1989), the difference is not surprising. Observed he- 
liocentric radial velocities for M stars (Maurice et 
al. 1987) in the same region as PSR J0045-7319 are 
similar to that of the pulsar. 

The results outlined here demonstrate that the 
B1V star identified by Kaspi et al. (1994) is the com- 
panion to PSR J0045-7319, and rule out a black-hole 
companion proposed by Lipunov et al. (1995). This 
makes PSR J0045-7319 the first dual-line binary ra- 
dio pulsar. The fact that no known radio pulsar has a 
black-hole companion means that we have no defini- 
tive proof that such systems exist, despite the exis- 
tence of objects such as Cygnus X-3 which might be 
expected to form a binary system comprised of a neu- 
tron star and a black-hole in the future (van Kerkwijk 
et al. 1992). This might be because binaries contain- 
ing black holes are particularly susceptible to coales- 
cence during the final common-envelope phase of their 
evolution. However, the expected number of observ- 
able pulsar black-hole binaries, given the statistics of 
the known pulsar population, is only of the order of 
one (Narayan, Piran & Shemi 1991, Lipunov et al. 
1994). A further implication of the companionship of 
the pulsar and B star is that the wind of the compan- 
ion is tenuous compared with those of typical Galactic 
B stars, as was suggested by Kaspi et al. (1994). Se- 
lection may play a role here, as higher stellar-wind 
densities might have rendered the pulsar invisible as 
in the case of PSR B1259-63 near periastron (John- 
ston et al. 1992). 
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